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, Abstract 

, ^ , ' We consider a 2 d.o.f. Hamiltonian system with one degree of freedom corre- 

sponding to fast motion and the other corresponding to slow motion. The ratio of 
! the time derivatives of slow and fast variables is of order < e <C 1. At frozen 

' values of the slow variables there is a separatrix on the phase plane of the fast vari- 

■ ables and there is a region in the phase space (the domain of separatrix crossings) 
where the projections of phase points onto the plane of the fast variables repeatedly 
cross the separatrix in the process of evolution of the slow variables. Under a cer- 

■ tain symmetry condition we prove existence of many (of order 1/e) stable periodic 
trajectories in the domain of the separatrix crossings. Each of these trajectories is 

1^ I surrounded by a stability island whose measure is estimated from below by a value 

' of order e. Thus, the total measure of the stability islands is estimated from below 

^ . by a value independent of e. The proof is based on an analysis of the asymptotic 

^ I formulas for the corresponding Poincare map. 

X: 

; 1 Introduction 

Many problems in the theory of charged particles' motion, the theory of propagation of 
short-wave excitations, and in celestial mechanics can be reduced to the analysis of 2 
d.o.f. Hamiltonian systems with fast and slow variables P,[2l[3]. One degree of freedom 
corresponds to fast variables, and the other corresponds to slow variables. The ratio of 
the time derivatives of slow and fast variables is of order < e <^ 1. To describe the 
dynamics in such systems one can use the adiabatic approximation constructed as follows. 

Consider the fast system, i.e. the system for the fast variables at frozen values of 
the slow variables. This is a 1 d.o.f. Hamiltonian system involving the slow variables as 
parameters. Assume that for a range of values of the slow variables there is a region filled 
with closed trajectories on the phase portrait of the fast system. Then one can introduce 
"action-angle" variables in the fast system The "action" variable of the fast system is 
an adiabatic invariant (i.e. an approximate first integral) of the complete system: its value 
oscillates with amplitude ~ e on time periods of order ~ 1/e along a phase trajectory. To 
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describe approximately the variation of the slow variables, one should average the rates 
of their variation over the "angle" variable of the fast motion. This approximation of the 
real motion is called the adiabatic approximation [5j. The obtained 1 d.o.f. Hamiltonian 
system for the slow variables, involving the "action" variable of fast motion as a parameter, 
is called the slow system. 

Assume that for some range of values of the slow variables the "action" variable phase 
trajectories of the slow system are closed. Then in this range the motion can be described 
with V.I. Arnold's theorem on perpetual adiabatic invariance [6]: the "action" value along 
a trajectory perpetually undergoes only small, of order e, oscillations (under a certain 
generality condition). The considered region of the phase space is filled, up to a small 
residual set, with invariant tori, which are 0(e)-close to the invariant tori of the system 
of the adiabatic approximation. 

We shall consider cases in which there are separatrices on the fast system's phase 
portrait (see Figured]), and there is a region in the phase space (the domain of separa- 
trix crossings) where the projections of phase points onto the plane of the fast variables 
repeatedly cross the separatrix in the process of evolution of the slow variables. In this 
case the above described construction of the adiabatic approximation needs some modifi- 
cation, and the above formulated assertions on the adiabatic approximation accuracy are 
not applicable. 

We shall assume that a certain symmetry condition is valid: the areas inside the 
separatrix loops in Figure [H are equal. 

The main result of the present paper is as follows. Under certain generality conditions 
in the domain of separatrix crossings on every energy level there exist many, of order 1/e, 
stable periodic trajectories of period ~ Each of these trajectories is surrounded by a 
stability island, and the measure of this island is estimated from below by a value of order 
e. Thus, the total measure of the stability islands is estimated from below by a quantity 
which is independent of e. A stability island is a domain on an energy level bounded by a 
two-dimensional invariant torus. A stability island contains a discrete family of invariant 
tori contractible to the periodic trajectory. Introduce a "modified action" equal to the 
"action" for points inside the separatrix loops in Figure [H and equal to one half of the 
"action" for the other points. This "modified action" is a perpetual adiabatic invariant 
of the motion inside a stability island: its value along a phase trajectory perpetually 
undergoes only oscillations with amplitude of order e. Therefore, the stability islands are 
also islands of perpetual adiabatic invariance. 

We note that the existence of stability islands with total measure that is not small with 
e in the domain of separatrix crossings is quite unexpected. Visually, in many problems, 
this domain looks like a region of dynamical chaos (see [HIS])- It was shown in |9j that a 
single stability island cannot have a measure larger than 0{e). 

The existence of stability islands with total measure not small with e in the domain 
of separatrix crossings was established in flUi [TT] in the case of a Hamiltonian system 
with one degree of freedom and the Hamiltonian function slowly periodically depending 
on time. Here we generalize this result to 2 d.o.f. systems. Like in pU] [TTj. the proofs are 
based on the study of asymptotic formulas for the corresponding Poincare map. In [TO|[TT] 
these formulas were constructed with the use of asymptotic expressions for the jump of the 
adiabatic invariant at a separatrix crossing in systems with one and a half d.o.f. [T2 | [T3 t [T^ 
and for the variation of the "angle" variable between separatrix crossings in such systems 
[T5] . Below we use analogous formulas for systems with two d.o.f. [ISl [IZj- To find 
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linearly stable periodic trajectories we look for linearly stable fixed points of the Poincare 
map. The results on Lyapunov stability of the periodic trajectories (the fixed points) and 
on the existence of invariant tori surrounding the periodic trajectories (invariant curves 
around the fixed points) are provided by the Kolmogorov-Arnold-Moser (KAM) theory. 
An example of application of the present theory to the motion of a charged particle in 
the parabolic model of magnetic field in the Earth magnetotail is given in [TS] . 



2 Formulation of the results 

Consider a two degrees of freedom Hamiltonian system with Hamiltonian H = 
H{p, q, y, x), where q, e~^x are coordinates, and p, y are canonically conjugated momenta, 
£ > is a small parameter, H G C°°. The corresponding equations of motion are: 

dH dH dH dH 

oq op ox oy 

The variables p, q are called fast variables, while y, x are called slow variables. The 
Hamiltonian system for p, q at {y, x) = const is called fast (or unperturbed) system. 

Assume that at all considered values of the slow variables the phase portrait of the fast 
system looks as in figure [H On this phase portrait, there exist a non-degenerate saddle 
point C and separatrices /i,/2 . These separatrices divide the (p, g)-plane into domains 
Gi = Gi{y, x), i = 1,2, 3. Denote the value of H at the saddle point as hg = hs{y, x), and 
introduce E = E{p, q,y,x) = H — hg. On the separatrices E = We suppose that E > Q 
in G3 and < in Gi and 6*2 • The areas of the domains Gi, G2 are 5*1, 5*2 respectively, 
and S'3 = 5*1 + S2] Si = Si{y,x) . We will make the following symmetry assumption: 

A) The areas of Gi,G2 are equal: Si{y,x) = S2{y,x) = S{y,x). 

Put O = Q{y,x) = {S,hs}, where {■, ■} is the Poisson bracket with respect to the 
variables {y, x) : {/, g} = f^g'y - f^g'^. 
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Figure 1: Phase portrait of the fast system. 

In the unperturbed system, one can introduce canonical "action-angle" variables (/, </?) 
separately in each Gi. The "action" variable / = I{p, q, y, x) in the unperturbed system 
involves y, x as parameters: I = A/ (27r), where A is the area on the phase portrait of the 
fast system surrounded by a phase trajectory passing through the point {p, q). For all p, q 
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on the Hamiltonian level line {p,q : H{p,q,y,x) = h} in the domain Gi the function / 
takes the same value Ii{h,y,x). 

Introduce the continuous function I{p,q,y,x) as follows: 

Hp,(l,y,x) = I{p,q,y,x), if {p,q) G ^1,2(1/, x); 

Hp^ ^) = Hp^ a;)/2, if {p, q) e G^iy, x); 

i{p, q, y, x) = S{y, x)/{27i), if {p, g) G /i U /a U G. 
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Figure 2: Schematic picture of motion on the slow plane. 

The function J is an adiabatic invariant of the perturbed system. Far from the sepa- 
ratrix of the unperturbed system, its value along a phase trajectory is preserved with an 
accuracy of order e on time intervals of order 1/e |6]. When a phase point crosses a narrow 
neighborhood of the separatrix, / varies by a value of order 0{e) for the majority of initial 
conditions [16]. The approximation / = const along a phase trajectory is called adiabatic. 
In this approximation the dynamics of the slow variables is described by a Hamiltonian 
system with Hamiltonian Ho{I,y,x), where Hq is the function H expressed in terms of 
/, y, X. Consider the dynamics on an energy level H = ho. On the phase plane of the slow 
variables {y,x), the separatrix is represented by a curve T{ho) = {y,x : hs{y,x) = ho} 
(see figure [2]). This curve is called the uncertainty curve 0. It divides the plane into two 
domains: -D3, corresponding to G^ on the fast plane, and the other, D12, corresponding to 
Gi and 6*2 • In the adiabatic approximation, a phase point on the slow plane moves along 
level lines I{ho,y,x) = const. As the area inside the separatrix loops varies, a phase 
trajectory can cross the separatrix. Fix u = 1 ot u = 2. Let B^{Io) be the adiabatic 
trajectory with / = Jq on the plane of the slow variables for phase points passing through 
Gj^: 

B,yiio) = {y,x: h{ho,y,x) = 2/o, {y,x) G ^3; I^{ho,y,x) = Jq, (y, x) G L'1,2} (2) 
Let us make the following assumption. 

B) There exists I^, such that the curve B,^{I^) is closed and intersects the uncertainty 
curve T{ho) at two points; at these points Q{x,y) 7^ 0. On By there are no stationary 
points of the Hamiltonian Ho- 

If assumption B) is valid both for u = 1 and v = 2 (with the same /=„), then in the 
adiabatic approximation, the projection of the motion with I = Iq onto the slow plane 
looks as follows. In the domain corresponding to G3, a phase point moves along the curve 
l3{h,y,x) = 2/*. After crossing T{ho), it follows one of the curves Iy{ho,y,x) = [in 
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general, these two curves are different]. Then the phase point crosses r(/io) again and 
continues its motion along l3{ho,y,x) = 21^. 

Assumption B) implies that there exists an interval H of values of /, such that for all 
Iq E D the curve B^{Iq) has the property B). In the exact system, multiple separatrix 
crossings occur in this interval of values of /. Computer simulations demonstrate diffusion 
of the adiabatic invariant and chaotic dynamics in this region, see, e.g., [TS]- 

The main result of the paper can be formulated as follows. 

Theorem. Under additional generality conditions, in the domain {p, q,y,x : / G S} 
on the energy level H = ho there exist more than /e stable periodic trajectories of 
system ([T]). Each of these trajectories is surrounded by an invariant torus bounding a 
domain of (3-dimensional) phase volume greater than C^^e. Inside each of these domains 
the variation of / is smaller than C^e. 

Corollary. The total measure of these domains is larger than C^^'^ = Cf ^6*2"^. 

Here Cj are positive constants independent of e. The generality conditions mentioned 
in the theorem are formulated below, in Sections and [71 

The domains bounded by invariant tori are domains (or islands) of perpetual adiabatic 
invariance. They are also called stability islands. Projected on the phase plane of the 
fast variables (p, q), the motion along the stable periodic trajectories under consideration 
looks as follows. A phase point starts moving in region G3, then crosses the separatrix and 
gets into region Gu, then crosses the separatrix again, comes back into G^, and returns 
to its starting point. 

The stable periodic trajectories mentioned in the theorem correspond to stable sta- 
tionary points of the return map generated by the original system. This map is written 
in the principal approximation in Section [51 

Remark. The result of the theorem is valid also if the phase portrait of the fast 
system is of the kind shown in figure [H but not for all x, y. It is enough to require that 
the system has such a phase portrait at the time of separatrix crossing calculated in the 
adiabatic approximation. 



3 Adiabatic and improved adiabatic approximations 

In the fast system, action-angle variables /, (y9mod27r are introduced separately in each of 
the Gi domains by a canonical transformation of variables. The corresponding generating 
function W{I,q,y,x) contains y,x as parameters [for brevity, we omit subscripts i]. We 
take this function in the form 

W{I,q,y,x)= f V{I,q',y,x)dq', (3) 

where V is the value of the p-variable along the trajectory with the prescribed value of 
the action /, and qo{I,y,x) defines a curve in Gi transversal to the phase trajectories. In 
the new variables the Hamiltonian has the form H = Hq{I, y, x). 

Now make a canonical transformation of variables {p, q, y, x) 1— > (/, (f, y, x) with the 
generating function ye~^x + W{I ,q,y,x). The canonically conjugated pairs of variables 
are (/, and {y,e~^x). Formulas for the transformation of variables are: 

i^ = dW/dI, p = dW/dq, x = x + edW/dy, y = y + edW/dx. (4) 
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In the new variables, the Hamiltonian H has the form 



H = Ho{I, y, x) + eH,{I, <f, y, x) + 0{e^), (5) 

where 

^ _dHdW _ dHo dW 
dy dx dx dy 

In the adiabatic approximation, the dynamics is described by the Hamiltonian Hq. In 
this approximation I = const along phase trajectories. 

One can also construct a canonical, close to the identity, transformation of variables 
{I,ip,y,x) ^ {J,ip,y,x) in order to make the terms of order e in the Hamiltonian inde- 
pendent of the phase (see, for example, [IZj). In the new variables, the Hamiltonian takes 
the form: 

H = Ho{J,y,x)+eni{J,y,x)+e'^n2{J,'ip,y,x,e), ni = {Hi), (7) 

where the brackets denote averaging with respect to if. 

In the improved adiabatic approximation, the dynamics is described by the Hamilto- 
nian Ho{J, Y, X) + £:7ii( J, Y, X). In this approximation J is an integral of motion. With 
an accuracy of order e^, the following formula for J is valid (see [16]): 
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2n Jq \2 

The integrals here are calculated along a phase trajectory of the fast system passing 
through the point {p,q); cr is the time of motion along this trajectory starting from this 
point, T is the period of motion. The function J is the improved adiabatic invariant. In 
the complete system far from separatrices its value along a phase trajectory is constant 
with an accuracy of order e"^ on time intervals of order e~^. In what follows, it is convenient 
to use the function J{p, q, y, x) equal to J in the domains Gi, G2 and equal to J/2 in the 
domain G^. 



J = J{p,q,y,x) = I + eu{p,q,y,x), (8) 



7^ r'?^6a]M- r^d.idt 



\dy Jo dx J Jq \ dx Jq dy 



+7^ r ('^-a]{E,h,}da. (9) 



4 Description of the separatrix crossing 

On the phase plane P^' of the slow variables {y, x), the separatrix is represented by a curve 
T{ho) (see figure [2]). Let By{I) G P'^' be a curve (see ([2])) corresponding to / such that 
/ e S. According to B), B^{I) is divided by T{ho) into two parts. Consider two segments 
1S3 and iSjy on P'^' separated by T{ho), such that each of the segments crosses transversely 
all the curves By{I) with / G H at a distance of order 1 from T{ho)- Points on the fast 
plane corresponding to 1S3 belong to G3, and the points corresponding to Sy belong to 
Gtj. Fix an interval Sq G S, assuming that the endpoints of Sq and S are different. 

Let the motion start at t = in a point M'^^\p^^\q'^'^\y'^'^\x'^'^'^)^ such that 
i/(p(o)^^(o)^ ^(0)^^(0)) _ ^Yie plane P^' {y^^\0^) G S^. Hence, (p(°),g(°)) G 
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G3(z/^°\a;(°)). Let / = /(°), J = J^^\ ^ = at this point; /(°) G Sq. In the adiabatic 
approximation the trajectory of the slow variables is B{I^^^), the motion is described by 
a Hamiltonian system with Hamiltonian Ho{I^^\y,x). Let and denote 

the points of intersection of B(I^^^) and r(/io); the sign "— " corresponds to passage from 
G3 to and the sign "+" corresponds to passage from Gi, to G3. Below r_ and r+ 
denote the slow time moments of separatrix crossing in this approximation. In the exact 
system, in the process of evolution the phase point on P**' approaches the curve T{ho), 
and accordingly, on the plane {p, q) it approaches the separatrix. We are interested in the 
dynamics of phase points that are being captured into G^, after the separatrix crossing. 
For such a phase point, the first crossing of the Cg-axis in G^ (see figure [I]) occurs near 
C at time r = r_ + 0(e). [We assume for simplicity that coordinate axes in figure [Hare 
the principal axes of the saddle point.] Denote the value of E at the point of this crossing 
as h^°\ We introduce t]^^) = 1 - \h^^yeQ^\ . Here e_ = e(y_,x_) > 0. 

After the separatrix crossing the phase point in the plane of slow variables {y, x) moves 
towards S^. When it crosses S^, the projection of the phase point onto the fast plane is 
deep inside region the G,^. Denote the value of J at this time moment as J^^\ 

Then the phase point starts again approaching the separatrix. At r = r+ + 0{e) the 
phase point crosses Cg-axis in G^ near the point C for the last time before entering G3. 
Denote the value of E at the point of this crossing as h^^\ We introduce r]^^^ = Ih^-^^eQ-f-] . 
Here 6+ = ©(?/+, x+) < 0. 

After crossing T{ho), the projection of the phase point onto the plane {y,x) crosses 
again the segment 1S3. Let J^'^\ ip^'^^ denote the values of J, ijj at this time moment. Then 
the phase point approaches the separatrix again, crosses it and gets captured into Gi, I = 1 
or / = 2. Let E = h^'^^ at the first crossing of the Cg-axis in Gi near C. This crossing 
occurs at time r = r_ + Tq + 0(e), where Tq is the slow time period of motion along the 
trajectory / = J'^°) in the adiabatic approximation. We introduce rf^"^^ = 1 — \h'^'^'> / eQ ^ \ . 
We are interested in the dynamics of the phase points for which / = i/. 

5 The return map 

In the energy level H = Hq the segment 1S3 is represented by a piece of a two-dimensional 
surface {p,q,y,x : H{p,q,y,x) = ho,{y,x) G S3}. This piece can be parametrized by 
variables J,^. The corresponding Poincare return map M : (J^^\ijj^^^) — > {J^'^\tjj^'^^) 
produced by trajectories that pass through G^ is symplectic. Its stable stationary points 
correspond to stable periodic orbits of the original problem, of period approximately equal 
to Tq/e. It is convenient to study these stationary points using a different set of variables, 
namely to consider the map M : (J^'^\ ri^^\ u) (J^'^\ri^'^\l). 
The map M is the composition of two maps: 

: (J(°),r/(°),i/) (J«,r/W,i/), M^^) : (J«,r/«,z/) ^ {J^'\r]^'\l) . 

The results of [161 [H] give the following formulas for M*^'^^ Suppose that 

ry(°),r^(^ry(^)6[crM-cr^ (10) 

Then 
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j(i) = j(o) _ (27r)-i£a_0_ln(2 sinTrr^^")) + {2n)-h{d^^^ - 4,^/2) + 

(27r)-i£0_(r7W - l/2)(6,,_ - 63,-/2) + 0{e^/Hne) , (11) 

^(1) = {r/W + £-i$i^^(J«,£) + 0(£i/3ln"i/3£)}, (12) 

j(2) _ J«-(27r)-^£a+0+ln(2sin7rr/W)-£(27r)-^(rf^,+ -d3,+/2) + 

(27r)-i50+(l/2 - r^(i))(6.,+ - 63,+/2) + 0{e'/'\ne) , ' (13) 

r^(2) = {r7«+£-i$2(J^'\£) + 0(£nn-i/='£)}, (14) 

l = u, if < {(£-^<l'2( J('\ e) + r/(^V2} < 1/2, (15) 

ly^u, if 1/2 < {(£-i$2( J('), e) + r/«)/2} < 1. (16) 



Here { ■ } denotes the fractional part, 

a± = aiy±, x±), bj^± = bj{y±, x±), dj^± = dj{y±, x±), 
j = 1, 2, 3, 63 = 61 + 62, 4 = + (^2- 

For a one has a = l/y/—g, where g is the Hessian of E at the point C. The values (ii^2 
are defined as the main terms of the expansions (cf. P^ ) 

27iu{p, q, y, x) = dj + 0(A/|£]'ln|i?|), 

where u is given by ([9]), and the point {p,q) is on the Cg-axis (figure [T]) in Gj near C. 
Values bj come from the expansion of the period of motion in G^, j = 1,2 near the 
separatrix: 

Tj = -a ln\E\ + bj + 0(Mn|^|). 

In (m, 

^[^\j,e) = ^j2^{ut\lY^{r),X^{r))+eu:t\lY^{r),X^{r)))dr, (17) 
^2{le) = - [ ' iJ^\2J,Y3iT),XsiT))+eJ^\2J,YsiT),XsiT)))dT. (18) 

Here uJq^ and cjj-''' are defined as uJq'' = dHo/dJ, tup'' = dHi/dJ, with Ho,Hi calculated 
in the region Gj, and {Yj,Xj) is a solution of the Hamiltonian system with Hamiltonian 
H = Hf){J, y, x)+eHi{J, y, x) on the energy level H = Hq. The values t'£^ are the slow time 
moments when the phase point corresponding to this solution arrives to the separatrix. 
At these moments SjiY^X) = 2ttJ. [One can take any of such solutions, they differ by a 
time shift which does not change the value of the integrals in (fT7|) . (ITSl) ]. Conditions ffTOl) 
are understood inductively: if and ri^^\ found according to (ITT]) . ffT2l) 

lies in segment [cj^\ 1 — c]"^], then (IT^ is valid, and analogously for r^^^^ and expressions 
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Equations (ITT]) . (fT3l) follow directly from the formula for the jump of the adiabatic 
invariant at a separatrix [16\ for the case when the areas of domains Gi and G2 are equal. 
Equation (fT2|) follows directly from the formula for phase change between separatrix 
crossings [TT]. Equations fll4l) - fll6l) follow from the result of [T7] in the following way. 

Let be the values of t and E when the phase point crosses the axis Cp near 

C for the first time after exit from G^,. At time t+ the phase point is in G3 on the 
positive or negative part of the axis Cp. Let h_ be the value of E when the phase point 
crosses the same part of Cp for the last time before exit from G3. Introduce notations: 
C = 2ll^' C = According to [17], 

( + ( = _L$2( + 0{e^/Hn-^/^e) mod 1, (19) 

provided (,( ^ [^2^)1 ~ '^2^]- The values h^^^ (see Section Hj) and /i+ are related as 
h+ = + e\Q+\ + 0(^3/2) [I6J. Therefore, ( = i(l - 77(1)) + 0{e^/^). The values 
(see Section Hj) and /;,_ are related as /i_ = /i*^^^ + 2eQ- + 0(£:'^/^) if Z = z/, and as 
h_ = + 50_ + 0(5^/^) if / 7^ I/. We have / = provided that 1/2 < ( < L Consider 
this case. We obtain that ( = i(l + r/(2)) _^ 0{e^/'^). From ([19]) we find 

^(1 - r^W) + i(l + 7^(2)) = 1$2(J(2)^ e) + 0{e^/Hn-^/h) mod 1. 

Therefore, 

r/(2) = r^(i) + 1$2( J('), e) + 0{e^/Hn-^^''e) mod 2. (20) 

Then condition 1/2 < ^ < 1 implies equation ( ITSil . 

Now consider the case I ^ u. In this case < C < 1/2. We get that C, and 7]^"^^ are 
related as C = + 0{e^''^). From ([I9]) we find 

1(1 _ ^(D) + 1^(2) = l$2(j(2)^ ^) + 0(£^/=^ln-i/3£) mod 1. 

Therefore, 

^(2) _ ^(1) ^ 1$2( j(2)^ e) - 1 + 0(ei/=^ln''/^e) mod 2. (21) 

e 

Combining (!20|) and (!2T|) we obtain equation (fT4l) . Condition < C < 1/2 implies (fT6|) . 



6 Linearly stable solutions 

6.1 Equations for stable stationary points 

Linearly stable periodic solutions of the original problem of period approximately equal 
to To/e correspond to linearly stable fixed points of the map M. These points are defined 
by the following set of conditions: 

J(2) = J(0)^ ^(2)=^{0)^ 1 = ^^ (22) 

|trM'|<2, (23) 
where M' is the linearization of M at the fixed point. 
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Put $1 = ^^^K Introduce also the following notations: 

jW=£^(^), 5 = 0,1,2; (24) 

a±0± = =F27ra± , di,,± - d3^±/2 = -2nd± , 6±(6^,± - 63,±/2) = 27r6±. 

Suppose that {^^^\ri^^^) = {^,1]) correspond to a stationary point of the map M. We 
rewrite equations ( l22l) using formulas (|TOl)-(|T^. notations (Ell), and neglecting terms 



rjW - r] + -f° aAn{2simrr]) + -f^d^ - liiv - ^)b- = £"^<l>i(£^, e) mod 1, (26) 



a_ln(2sin7rr/)+rf_-(r7-l/2)6_ = a+ln(2sin7rr7(^))+rf++(l/2-r/(^))6+, (25) 

- ry(i) = e-i$2(£^,e)modl, (27) 

£-^<l>2(e^,£) e (0, l)mod2, (28) 

r7,r7« G (q\l-cr'). (29) 



V 



(1) 



The stability condition | trM' | < 2 can be written using (fTT|) - (fT4l) [we again neglect 
terms 0{e^/Hn~^/^e)]: 

-4<g<0, (30) 

Q={U2~ «i)(72 + 7?) - U2uahl (31) 
Ml = a-TT cot TTi] — b_, U2 = a+VT cot vrr]'-"'^'' — b+. 

The exact set of conditions for a stable periodic solution differs from ( I25l) -( l29l) by terms 
0{e^/Hn-^/^e) that can be differentiated in ,1] , t]'- > without changing their smallness 
order. Below, we find non-degenerate solutions of system (!25|) - fl29l) . According to the 
implicit function theorem, for small enough e, each of these solutions corresponds to a 
solution of the exact system. 

6.2 Existence of stable stationary points 

We assume that the following generality conditions are valid. 
HI. In Hq there exists an interval of values of I such that 

7?(/) ■ 72°(/) ^ , ^(7?(/)/72°(/)) ^ . (32) 

In what follows, we consider only those values of eC, that belong to the interval defined by 
HI. 

Equation (!25|) has the form 

F{v; b-, dJ) = F(r/(^); a+, 6+, d+) 

with F{r]; a, b, d) = a\n{2 sinvr?]) + d — b{r] — 1/2), t] G (0, 1), a > 0. We have F — oo 
as ?7 — > 0^ or 7] — > 1~. Moreover, F has a unique maximum F^ax = m{a,b,d) on the 
interval < t] < 1. We put m± = m{a±, b±, d±). 
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Figure 3: Set of points {r],!]^^^) satisfying fl25l) . a) m_ > m+; b) m_ = m+; c) m_ < m+. 



The set of points (t],?]^^-*) G (0,1) x (0,1) satisfying (!25|) is shown in figure [31 Con- 
sider for definiteness the case m_ > m_|_, 7^72 < 0. The other cases can be studied 
analogously. 

Let in figure [3] A and B have coordinates (0, 0) and (0, 1) respectively. Consider 
condition fl30|) on the curve AB. Near the point A both ui and U2 are large positive, 
and therefore Q is large positive. Near the point B both ui and —U2 are large positive, 
and therefore —Q is large positive. Therefore, on the segment AB there is an interval 
L = L{e^) (or several intervals) where condition (130|) is valid. If ci is sufficiently large, 
the relations ( l29i) are also valid on L{eC,). 

Now we look for solutions of the system 



(r7,r7W)eL(£0, 

?7 + 7°a_ln(2sin7rr7) + 7?ci- - -fib-{r] - 1/2) = £:"^<l)i(£:^, £:) modi, 

- r]^'^) = s'^^^(^s^,e) mod 2. 



(33) 
(34) 
(35) 



Relations fl33|) - fl35|) are equivalent to fl25l) - fl3T|) . Relation fl28l) holds because, according 

= 1] mod 2 and t] G (0, 1) . Validity of the other relations 
follows from the definition of L{e^) and equations flSlj) . ([35D- 
System fl5^ - fl5S]) can be interpreted geometrically as follows. Consider the curve 
A(e^) C = {(si mod 1, S2 mod 2)}, which is the image of L{eC,) under the mapping 

{ri,r]^^^) (si,S2) = {r]'-^^ - + 7?a_ln(2 sinvr?]) + 7?rf- - liiv - 1/2)^-, rj - rj^^^) 

defined by the left hand sides of (IMjl . (1351) . Consider also the point 

e-^$(e^,e) = (£"^<l>i(£^,£)modl, £^^$2(6^, e) mod2) 

on T^. As ^ varies, the curve A(£:^) is moving slowly on the torus. At the same time, 
the point is moving fast, with its velocity vector 7(5^, e) = (71(5^, e) , 72(5^, e)) 

varying slowly. Solutions of system (!33l) - (!35l) correspond to values of "time" ^ when point 
e~^(^{e^,e) crosses the curve A(e^). 

It follows from assumption HI, that there exists 1^ such that the ratio 7i(-^c)/72(-^c) 
is irrational. Then in the interval ,^ G {Ic — C3 ^/^; Ic + /^) on every segment of 
length C4 of "time" ^ the point e~^^{eE,iS) crosses transversely the moving curve A(£:^) 
at least once. Indeed, the point's trajectory is close to a dense winding of the torus with 
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irrational angular coefficient 7i(-^c)/72(-^c)j and the curve A(Jc) is transversal to the vector 
7°(/c) = (7i (-^c), 72(-^c)) at least at one point. [The latter is because A(Jc) is not straight.] 
Therefore, there exist more than c^^c^^ je — 1 > C^^ /e intersection points. Each 
intersection point corresponds to a linearly stable stationary point of the return map. In 
Section [S] we obtain an asymptotic expression for the number of intersection points. 

7 Stability islands 

Let ^r, rjr, T]^^ be one of the solutions of system (1331) -(l35l). found in the previous section. 
One can rewrite the map M in variables ^ = ^ — Vi expand its right hand sides in e 
supposing f ~ 1, and neglect terms 0(£:^/^ln~^/^£:). Thus, one obtains a map that does 
not contain e and depends on the parameter = eE,r- This map has a stationary linearly 
stable point ^ = 0, = rjj.. According to Kolmogorov-Arnold-Moser (KAM) theory, 
under certain conditions [absence of resonances up to the 4-th order, non-zero coefficient 
in the normal form] the stationary point of the complete map is surrounded with a Cantor 
family of invariant curves (see, e.g. [5|). These latter conditions of KAM theory in the 
case under consideration can be written as one condition, namely /(/r, Tjr) 7^ 0, where / 
is a continuous function. We assume that the following condition is valid. 

H2. There exist J*, r/*, r/l^^ (r/*, r/l^^) G L{h) such that /(/*, r]^) 7^ . 

Then for / close to |/ — /*| < Cg \ there is a segment of the curve L{I) where 
|/(J, Tj, r/^^-*)! > Cg ^ Without loss of generality one can assume that this segment coincides 
with L[I) [otherwise, we can shorten L[I)]. According to KAM theory, for — /*| < c^^ 
there is an invariant curve at distance bigger than Cg ^ from the point ^ = 0, r] = rjr. We 
call the domain surrounded by this curve a stability island of the Poincare map. The area 
of this island is larger than cf^. Returning from ^, i] to the original variables ^, 77, we find 
that the area of the stability island is estimated from below as Cj^e. 

The boundary of the stability island (an invariant curve) is a section of an invariant 
torus on the energy level of the original system. The volume of the domain surrounded 
by this torus is larger than C^^e. The variation of / inside this domain is smaller than 



8 Distribution of periodic orbits 

Let N(Ic, A,e) be the number of stable periodic orbits under consideration with /^"^ G 
[Jc, /c + A] C Sq. According to Section [6^ this number is equal to the number of "time" 
points ^ e [Ic/£, (-^c+A)/^] when the point e~^^{e^, e) crosses the curve A{e^) on the torus 
= {(si mod 1, S2 mod2)}. For values of Ic such that 7i(/c)/72(-^c) is irrational one can 
calculate p(/c) = limA^o lime-*o(^A^(^0) -^c)/A). We give the result omitting the standard 
argument based on ergodicity of winding of a torus with frequencies (7°(/c), 72(-^c)) (cf. 
|19j). Let / be a natural parameter on A(Jc), and n{l, Ic) be a unit normal vector on A(/c). 
Introduce the notation 



where (7°, n) is the Euclidean scalar product of 7° and n. Hence, x is an "absolute" phase 
flux of the vector field 7° across the curve A: for the phase flux across any element of A 




(36) 
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we take its absolute value. The value of x(/c) equals the product of the length of vector 
7^^ and the full length of the projection of the curve A(/c) onto the normal to vector 7°. 
Then 

p{Ic) = \x{Ic) (37) 

(the coefficient 1/2 in this formula appears because the surface area of the torus equals 
2). 

Let iV(So, e) be the number of stable periodic orbits under consideration with 7^°) e Sq. 
If 7id7i/d/ — 72d72/d/ vanishes only on a set of measure on Sq, then 

hmeN{Do,e) = l f x(/)d/. (38) 
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